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ABSTRACT 

We continue the study of time-dependent Hamiltonians with an isolated singularity in 
their time dependence, describing propagation on singular space-times. In previous work, 
two of us have proposed a "minimal subtraction" prescription for the simplest class of 
such systems, involving Hamiltonians with only one singular term. On the other hand, 
Hamiltonians corresponding to geometrical resolutions of space-time tend to involve mul- 
tiple operator structures (multiple types of dependence on the canonical variables) in an 
essential way. 

We consider some of the general properties of such (near-) singular Hamiltonian systems, 
and further specialize to the case of a free scalar field on a two-parameter generalization 
of the nuU-brane space-time. We find that the singular limit of free scalar field evolution 
exists for a discrete subset of the possible values of the two parameters. The coordinates we 
introduce reveal a peculiar reflection property of scalar field propagation on the generalized 
(as well as the original) nuU-brane. We further present a simple family of pp-wave geometries 
whose singular limit is a light-like hyperplane (discontinuously) reflecting the positions of 
particles as they pass through it. 



1 Aspirant FWO 



1 Introduction 



Defining dynamical transitions through space-time singularities entails a very large amount 
of ambiguity, both technically and conceptually, and, clearly, a deeper and more systematic 
understanding of gravitational physics is needed in order to address the issue with full 
legitimacy. Nevertheless, even in the absence of such understanding, it appears desirable to 
explore the range of possibilities presented by the problem of evolution across singularities. 

In several string theory approaches, holography maps the study of cosmological singu- 
larities to the study of quantum mechanics or quantum field theory with certain singular 
features: couplings may develop a singularity as a function of time ([I] and related work), 
the quantum field theory may live on a singular spacetime (p] and related work) or it 
may have a potential unbounded below [2]. In particular, the model of [1] is described by a 
quantum field theory of matrices on (the future cone of) the compactified Milne space-time. 
A similar model studied in [3] involves quantum field theory of matrices on the nullbrane 
space-time |1] and its singular limit, the parabolic orbifold The models presented in 
[6] are quantum mechanics models of matrices, making them simpler from some points of 
view. An important question is what happens in these field theories when the space-times 
they live on develop a singularity. While certain subtle questions related to the large N 
limit (A^ being the size of the matrices) have not been fully addressed yet, these models 
clearly motivate the study of field theory with (near-)singular Hamiltonians. 

The parabolic orbifold [7] , the nullbrane [8] as well as the compactified Milne space-time 
have also been studied as (a part of) backgrounds of gravitational theories, in particular 
as (a part of) string theory space-times. Most attempts to study the singular space-times 
among those just mentioned in string perturbation theory have failed because both the 
parabolic orbifold and the Milne orbifold exhibit divergences signaling large gravitational 
backreaction [3 El [lOl [H]- See, however, [12] for a closely related model in which the 
singularity is replaced by a phase with a condensed winding tachyon within perturbative 
string theory. Note that in the models we shall discuss, these singular space-times host 
a holographically dual quantum field theory, not a gravitational theory. In fact, in the 
present paper we limit ourselves to the study of free field propagation on a fixed (possibly 
singular) space-time background. Interesting questions related to possible non-gravitational 
backreaction are postponed to future work. 

In a previous publication [13] , two of us noted that evolution across space-like or light- 
like singularities appears to be often described by quantum time-dependent Hamiltonians 
with an isolated singularity in their time dependence. We then exposed the most conser- 
vative way to define a unitary quantum evolution corresponding to such Hamiltonians by 
modifying the singular time dependences to become distributions while keeping the opera- 
tor structure of the Hamiltonian unchanged. This approach is relevant when the transition 
through the singularity is dominated by a single term (single operator structure) in the 
Hamiltonian, and one can think of this way to define the transition through the singularity 
as a sort of "minimal subtraction". In the absence of further physical specifications, this 
procedure appears to be the most natural way to define evolution across singularities. 

However, in many geometrical contexts, another approach appears to be more natural. 
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Namely, one may want to resolve the singular geometry into a smooth space, and then try 
to take the singular limit in such a way that the dynamical evolution remains well-defined. 
This is non-trivial, since naive resolutions of a singular space-time will generically not lead 
to a well-defined dynamics in the singular limit. As we shall argue shortly, constructing 
such geometrical resolutions will typically take us outside the formal scope of [13], and it 
is the aim of our present paper to investigate what kind of mathematical structures the 
geometrical resolutions of space-like (light-like) singularities involve: in general, as well as 
in a few specific examples, which will be at the focus of our attention. 

In [13], we divised a kind of "minimal subtraction" prescription and described its par- 
ticular implementation for the case of a free scalar field propagating on the compactified 
Milne universe. Here, we would like to argue that the approach of [T3| does not lend it- 
self to a geometrical interpretation, and, therefore, should one be interested in geometrical 
resolutions of space-time singularities, a more general framework is required. 

To recapitulate briefiy, the metric of the compactified Milne universe is 

ds'^ = -dt^ + t^dx^, X r^x + 2T[, (1) 
and the corresponding free scalar field Hamiltonian is 

H=l^Jdx{nl + 0'^) dxc^\ (2) 

With this form of the Hamiltonian, the Schrodinger equation cannot be integrated through 
t = on account of the singularity of 

The idea of the "minimal subtraction" scheme of [13] is to keep the operator structure of 
the Hamiltonian unchanged and to modify the singular time dependences in (|2]) locally at 
t = by subtracting terms proportional to (possibly) resolved 5-functions and its derivatives 
such that the time dependences become well-defined in the sense of distribution^. Then, 
the Schrodinger equation can be integrated. Put differently, one can replace in ([2]) 
by its regulated version fii\t\(t,e) (with e being a regularization parameter), in such a way 
that, as e is taken to 0, fi/\t\{t,e) converges to a distribution Ti/\t\{t), and this distribution 
J-'il\t\{t) equals l/|t| everywhere away from t = 0. A possible choice is 

hmit^e) = -^ + 21n(/.e)-^^^, (3) 

with /i an arbitrary mass scale. 

With this approach, one obtaines a regularized version of the Hamiltonian ([2]), namely 

1 „ . f , ^,2,-/2 



H = -h/\t\{t,e) I dx {7r', + r)+--- (4) 



^This procedure bears a strong formal resemblance to the conventional renormalization of local field 
theories by subtracting local counter-terms, and it can be thought of (see [13j for further discussion) as 
renormalizing the time dependence of 
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e=0.05, fi=3 



Figure 1: Negative contribution around t = in regularised fi/\t\{t,e) ([3]) 

such that, as e is taken to 0, the evolution away from t = becomes identical to that arising 
from ([2]), and, furthermore, the system displays a well-defined (unitary) transition through 
t = 0. 

The "minimal subtraction" procedure we have just briefly re-stated, is a consistent evo- 
lution prescription in itself. However, a direct inspection of (jlj) shows that the regularized 
version of our dynamics does not admit a geometrical interpretation (nor should one think 
of its singular limit, albeit well-defined, as being geometrical). 

The problem with constructing a geometrical interpretation of (jl]) is that, since fi/\t\ (t, e) 
has an e —>■ limit as a distribution, the e —>■ limit of 

to 

J dthm{t,e) (5) 

-to 

must exist. Furthermore, as stated above, the e — limit of fii\t\{t,e) must equal l/|t| 
everywhere away from t = 0. For that reason, in order for the limit of ([5]) to exist, fij\t\{t, e) 
should be very large and negative somewhere in the e-neighborhood of t = so that the 
positive divergence from integrating is compensated in ([5]). This is clearly apparent 
in figure [H However, the coefficient of the kinetic term in the Hamiltonian of a field in a 
geometrical background comes from the square root of the determinant of the metric (and 
the coefficients of the inverse metric), and it needs to be positive (as is the function 
appearing in ([2])). 

For that reason, there appears to be a conflict between the demands of positivity for 
certain coefficients in the Hamiltonian arising if one pursues a geometrical interpretation, 
and negative contributions introduced by our "minimal subtraction" recipe. If one is to 
construct a geometrical resolution of dynamics on a singular space-time background, one 
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generally needs to relax the specifications of the "minimal subtraction" approach, and 
permit modifications in the operator structure of the Hamiltonian, as well as its time 
dependence, in the vicinity of the singular region. One will then typically end up with a 
situation where a few different operator structures in the Hamiltonian essentially contribute 
to the transition to the singular region: 



where Hi are time-independent operators and fi are time-dependent number-valued func- 
tions, e is a regularization parameter, and the implication is that, as e is taken to 0, some 
of the /j's may develop isolated singularities at a certain value of t, which we shall choose 
to be t = 0. It is the commutation properties of those different terms in the Hamiltonian 
that are responsible for divergence cancellation (rather than explicit negative contributions 
introduced through the "minimal subtraction" scheme of p3])- 

In what follows, we shall review some relevant properties of quantum time- dependent 
Hamiltonians involving multiple operator structures, and proceed with applying this range 
of techniques to constructing singular limits of dynamics in a few different geometrical 
backgrounds related to the null-brane space-time. In section 2, we start by discussing the 
evolution properties for Hamiltonians of the form ([6]). In section 3, we examine as an 
example the null-brane, which is the geometric resolution (through a parameter R, which 
plays the role of e in this case) of the parabolic orbifold. We introduce a generalized null- 
brane metric depending on continuous parameters {R,a,l3). The propagation of a free 
scalar field on this space-time will be described by computing its mode functions, derived 
by solving the wave equation exactly using WKB methods. Then we discuss properties 
of the mode functions of a free scalar field in the singular limit R ^ 0. We shall find 
that the limit of the mode functions exists for certain discrete values of the parameters a 
and p. We give a qualitative discussion of the results in subsection 3.3. In section 4, we 
introduce the light-like reflector plane, a simple yet curious space-time that can be used 
as a toy model for light-like space-time singularities. Section 5 contains our conclusions. 
Appendix A contains a discussion of the minimal subtraction prescription of |13] applied 
to the parabolic orbifold. In appendix B, we give some technical details related to focusing 
properties and Maslov phases. 

2 Time-dependent Hamiltonians involving multiple operator struc- 



In preparation for our analysis of geometrical resolutions and their singular limits, we shall 
review the quantum dynamics described by Hamiltonians of the form (jH]). Our ultimate 
question will be whether the e —>■ limit of the evolution operator corresponding to (jH]) 
exists. 

It is in general impossible to solve the Schrodinger equation corresponding to the Hamil- 
tonian ([6]). The familiar symbolic solution for the evolution operator U{ti,t2) involves the 




(6) 



tures 
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time-ordering symbol T: 



U{ti,t2 



t2 

-i j dtH{t) 



(7) 



The above representation can be further transformed in an instructive way using a 
technique known as the Magnus expansion [T3]. The operator U belongs to the group of 
unitary operators on the Hilbert space, and the Magnus expansion can be thought of as 
an analog of the Baker- Campbell-Hausdorff formula for finite-dimensional Lie groups (the 
latter is discussed in many textbooks on group theory, for example, in [15]). The expansion 
can be symbolically written as: 



exp 



t2 



-i / dtH{t) + r]i / dtdt'[H{t),H{t')] 



+ir]2 / dtdt'dt"[H{t),[H{t'),H{t")]] 



with some numerical coefficients t]i, f]2, ■ ■ ■ (their values will not be important for us, and 
it appears they can only be derived recursively [H]). The key property of the above ex- 
pression is that the higher order terms are entirely expressed through higher order nested 
commutators of H{t) at different moments of time. 

Even though, in a completely general setting, the Magnus expansion is hopelessly in- 
tractable, it displays the broad range of opportunities for divergence cancellation in a 
singular limit of the dynamics described by ([6]). Namely, for the case of (Q, the Magnus 
expansion ([8]) will contain all kinds of combinations of the fi and their products, in such a 
way that, even if fi develop very strong singularities as e is taken to 0, the limit of f/(ti, ^2) 
may still exist. For example, even if all fi are positive, cancellations may still occur on 
account of the commutation properties of Hi. 

Should such cancellations take place, one may think of the e ^ limit of ([6]) as an 
operator- valued generalization of conventional distributions: just as ordinary distributions 
may contain singularities in a way that permits evaluating ordinary integrals, the Hamil- 
tonian will contain singularities in a way that permits evaluating the time-ordered 
exponential integral in ([7]). A systematic exploration of such generalized operator- valued 
"distributions" may be interesting to pursue, but lies outside of the scope of the present 
paper. 

There is a special case when the above analysis can be taken significantly further. 
Namely, it may turn out that, for all moments of time, the operator f/ of ([7]) belongs to a 
finite-dimensional subgroup of the unitary group of the Hilbert space. This situation has 
been described as a presence of a dynamical group (see jMl [T7] and references therein). For 
the Hamiltonians of the form ([6]), there will exist a finite-dimensional dynamical group if 
the set of nested commutators of ifj's closes on a finite-dimensional linear space of operators 
(which would serve as the Lie algebra of the dynamical group). Should that happen, one 
would be able to use the the closed resummed version of the Baker-Campbell-Hausdorff 



5 



formula for finite-dimensional Lie groups (see, for example, [15]) to treat the Magnus ex- 
pansion, or, alternatively, the Schrodinger equation can be reduced to a finite number of 
ordinary differential equation^ describing the evolution on the finite-dimensional dynamical 
group manifold [Ml HJl US] ■ In practical terms, one can choose a particular low- dimensional 
faithful linear representation of the dynamical group furnished by matrices M, and write 
down the Schrodinger equation in this representation: 



.dM{t,to) 
' di 



^{H{t))M, M(to,to) 



(9) 



where is a homomorphism from Hilbert space operators onto the representation furnished 
by M. (This is a finite-dimensional system of ordinary differential equations.) Given the 
solution for M{t,to), one can reconstruct the original evolution operator as ip~^{M{t,to)). 

While the specific examples of quantum dynamics discussed in this paper will be con- 
structed using a somewhat unconventional application of WKB methods, the fact this 
"double-semiclassical" analysis is possible reflects an underlying dynamical group structure 
inherent to the systems we shall work with. After explicitly computing the mode functions 
encoding the dynamics, we shall investigate (in this greatly simplified setting) the existence 
of singular limits. It will then be possible to circumnavigate the formal complications intro- 
duced by the non-commuting structures in ([6]), and examine the limiting case of evolution 
on a singular space-time background. 

3 Generalized null-brane space-times and singular limits 

In this section, we study quantum dynamics of a free scalar field propagating on the 
parabolic orbifold and some geometrical resolutions thereof, including the null-brane. 



3.1 The null-brane geometry, the parabolic orbifold and generalizations 

The null-brane spacetime, originally introduced in [4j, was studied in the context of per- 
turbative string theory in [SJ [231 121]; a matrix theory description was provided in [3]. 

Consider Minkowski space-time in light-cone coordinates ds"^ = —2dx^dx~ + dx"^ + dz"^. 
The null-brane is obtained by identifying 



X ' 
X 

x~ 



exp{2nnj') 



X ' 

X 

x~ 



z + 2'nnR . 



n G Z . 



J 





1 
1 



• (10) 



In the -R ^ limit the null-brane reduces to the parabolic orbifold [51, ^7] times the line 
labeled by z. In this sense, the null-brane is a geometrical regularization of the parabolic 



■^The analytic power of the dynamical group approach does not appear completely clear or fully explored. 
It certainly does apply to all linear quantum systems; however, in that case, the conventional WKB analysis 
would suffice. Beyond linear systems, the relevant finite-dimensional subalgebras of Hermitean operators 
may be difficult to construct and/or classify. Nevertheless, some non-trivial examples of dynamical groups 
for quantum- mechanical systems do exist (see, for example, [I9j). 
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orbifold. In Rosen coordinates, the metric on the parabohc orbifold is given by 



ds^ = -2dy+dy- + {y+fdy'^. (11) 

In [8], the null-brane geometry was discussed in two coordinate systems, 

ds^ = -2dy+dy- + du^ + {R^ + {y+y)dy^ + 2Rdydu; (12) 
ds^ = -2dx^dx- + dx'^ + {R^ + x'^)de^ + 2{S:^dx-xdS:+)de, (13) 

related to Minkowski coordinates by 

x+ = ?/+, x = yy^, X' =y' + ^V^V^, u = z - Ry; (14) 

_ 1 _ z 

x^ = x'^, X = X + 9x'^, x~ = x^ + 6x -\ — O'^x'^, 9 = —. (15) 

2 R 

Unfortunately, neither coordinate system is fully satisfactory for studying the R —y limit 
of dynamics on the null-brane0 On the one hand, the y-coordinates are not globally defined 
since they are singular at = for any R. On the other hand, the x-coordinates, which 
are nonsingular for R ^ 0, do not have an i? — limit even away from the parabolic 
orbifold singularity, as the determinant of the metric is —R^ everywhere. 

Therefore, we now introduce new coordinates that interpolate between the x-coordinates 
(for small x"*") and the ^/-coordinates (for large y~^): 

X+ = x+ = ?/+; (16) 

. 1 i A R^ \ R^ y+u^ 

X- = X---— l-——-—-^)=y- + 



X 



2x+ \ ^R2 + ^x+yYJ 2 {R^ + {y- 

Rx y^u 



\2 ' 



n ^ f. R^ \ Ru 



X+ \ i?2 + (£+)2 J ^ ' i?2 + (y+)2 

In this coordinate system, the metric has determinant —{R^ + (X+)^). 

The metric of the null-brane space-time, written in our new coordinates, can be naturally 
generalized to a two-parameter family of metrics, which we will call generalized null-brane. 
The family is labeled by parameters a and /3, the original null-brane corresponding to 

a = 3,/3 = 2: 

ds"" = ~ ""l (dX+? -2dX+dX-+ , ^^^^ dX-^dO+jR^ + (X+)2) d(d'+dX\ 

(17) 



Following the discussion in the previous section, with R playing the role of e, we would like to phrase 
the dynamics in terms of a Hamiltonian that has the structure H — fi{t, R)Hi in which fi(t, R) is 
regular in t for R ^ and regular away from t = for R = 0. The terms that appear in the Hamiltonian 
can be easily deduced from the inverse metric. 
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3.2 Free scalar field on the generalized null-brane 
3.2.1 General dynamical preliminaries 

We consider a free scalar field on the generalized null-brane metric. After Fourier trans- 
forming with respect to X~ and G, 



fee 



'It 



and suppressing the indices k- and ke, the action reads 



ke 

i2 



ik_ 
~2 



d 



ftA aX^R^k^ kek_pXR 



2 2(i?2+(X+)2) 2(i?2 + (X+)2)2 (i?2 + (X+)2)3/2 



(.19) 



Denoting dx+(p as 0, the wave equation reads 



iX- 



d 



+ 



pXRke 



2(i?2 + (X+)2)- 2k. ' (i?2 + (x+)2)3/2 



(20) 



^2 



+ 



+ 



m 



2fc_(i?2 + (X+)2)- 2 (i?2 + (j^+)2)2- 2A;_^ 

As the Lagrangian is first order in time derivatives, we have to deal, in principle, with con- 
straints when deriving the corresponding Hamiltonian. One can take a shortcut interpreting 
as a canonical coordinate and tt = ik-^^/R'^ + (X+)20* as its conjugate momentum. The 
Hamiltonian reads 



H 



EE 

fce fc-T^O 



dX 



TT 



2 (i?2 + (X 



A-2 



fm^ i ^2 ipXRke 
2^ + 2^ ~ (i?2 + (x+)2) 



3/2 



2fc_ i?2 + (X+)2 2 (i?2+(X+)2)' 



(21) 



It is manifestly of the form H = fi{t, R)Hi, in other words, it belongs to the class of 
Hamiltonians we single out in section 2. 

We now show that the Hamiltonian ( l2Ti) leads to a finite-dimensional dynamical group 
structure of the type discussed in section 2. The canonical variables 7r(X, t) and (f){X, t) ap- 
pear in four combinations J vr d\(j)^ J 7r0, J 7rX0 and J irX'^cf). The commutation relations 
for such operators are given by 



TT A{X) 4>dX, 71 B{X) (f) dX 



n[A,B](^dX, 



(22) 
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reducing the commutator algebra to that of {dx, 1, X, X^}, which closes after the addition 
of {Xdx, dx}- Equivalently, one can form the standard (single degree of freedom) creation 
and annihilation operators a and out of X and dx to conclude that the Lie algebra of 
the dynamical group is spanned hj n = a^a + 1/2, a^^, a^, a\ a and /. Given the inclusion of 
powers of the creation operator up to a^^, it is not surprising that the corresponding algebra 
has become known as the two-photon algebra, or h^, and has been featured in discussions 
of quantum optics, and squeezed states in particular (see, for example [IS]). A complete 
formal analysis of quantum dynamics on the two-photon group has been given in |18j . 

Following the general picture presented in section 2, we could use the two-photon group 
considerations of |r8j to reduce the question of free scalar field dynamics on the generalized 
nuU-brane to ordinary differential equations. In our present setting, however, one can 
perform these operations in a considerably more familiar guise. Namely, since the free 
scalar field is linear, solving for its quantum dynamics amounts to constructing a complete 
set of solutions to the classical wave equation. Furthermore, the classical wave equation 
turns out to be equivalent to the Schrodinger equation for a linear auxiliary one-dimensional 
quantum system. Because the auxiliary system is linear, its Schrodinger equation (i.e., 
the wave equation of the original scalar field) can be solved exactly by WKB methods. 
The latter effectively reduce the problem to ordinary differential equations (the classical 
equations of motion of the auxiliary linear system). Thus, one attains the same level of 
simplification as one would through performing the analysis of [I8l. One can refer to the 
above procedure as "double-semiclassical" analysis (there is an (exact) WKB procedure 
leading from a free quantum scalar field to the wave equation for the mode functions, 
and an (exact) WKB procedure leading from the wave equation for the mode functions to 
a one-dimensional auxiliary classical system). Note that both the "double-semiclassical" 
approach and the general dynamical group approach of [18j (which are essentially one and 
the same thing) are made possible by the fact that the metric of the generalized nuU-brane 
is a quadratic polynomial in the X-variable. 

For linear quantum systems, it is most common to work in the Heisenberg picture, 
instead of (equivalently) deriving WKB wave functions in the Schrodinger picture. One 
obtains the solution for the Heisenberg field operator as an expansion in terms of a complete 
set of mode functions Uy{X, X~^) (with y being a generic basis label) satisfying the classical 
equations of motion: 



If one demands the standard commutation relations for the creation-annihilation operators 
and a, the canonical commutation relation between vr and determine the normalisation 
of the mode functions (this is the analog for first order systems of the Klein-Gordon norm): 




(23) 



The corresponding conjugate momentum is 



7rk_,ke=^k.^R^ + {X+f [ dy u;{X, X+)a^{y) 



(24) 



S{X ~ X) = k_^J + {X+f j dy ul{X,X+)uy{X,X+). 



(25) 



9 



3.2.2 "Double-semiclassical" solution of the wave equation 

Denoting by t, fl20|) takes the form of an auxiliary Schrodinger equation with Hamiho- 



nian 



it P2 fjXRke kl a X^R^k. 



(up to a sign difference in the left hand side of (!20|) ). As the corresponding Hamiltonian 
( l26l) is quadratic in X, (120!) can be solved exactly by WKB methods. The starting point is 
the observation that the ansatz 

<f){X^,h\X2,t2) = A{h,t2)exp{-tS,i[Xi,ti\X2,t2]) (27) 
solves (I20D (with t) if 



Sci= I ' dt [px-n 

ti 



X=X,i{Xi,ti\X2,t2) 



{2i 



-2._^^^(t„t)^!^l^^. (29) 

Here, is the classical action with boundary conditions X{ti) = Xi, X{t2) = X2. More 
general solutions to (12UI) are obtained by integrating (17r|) over Xi, weighted by an arbitrary 
smooth wavepacket. 

A subtlety arises in this ansatz when the dynamical evolution reaches a focal point 
t2 = t*, where the classical action diverges, unless a certain relation between Xi ("the 
source") and X2 ("the image") is met. At such focal points, the differential equation for 
A{ti,t2) becomes singular. In that case, one solves the WKB equations away from t* 
and connects the solution by a phase jump at the focal point. The phase jump should 
be chosen such that convolutions of (p71) with a smooth wavepacket are continuous across 
the focal point. The general guidelines for this procedure are best familiar in the context 
of caustic submanifolds in geometrical optics (see, for example, [201 EI]) and the above- 
mentioned correction pre-factors have become known as the Maslov phases. We shall give 
some further details in appendix [Bl 

In order to compute the classical action, we first consider the classical equation of 
motion: 

X + a- -X = -^ — , 30a 

(i?2+t2)2 k. (i?2+i2)3/2 

X(ti)=Xi, X{t2)=X2. (30b) 

This equation is actually exactly solvable, and it has become known as the equation for 
"bending of a double- walled compressed bar with a parabolic cross-section" [22]. It can 
be reduced to a driven harmonic oscillator with constant frequency via substitution X = 
A/i?2 -)- t'^xivi't)): taking 77 = arctan(t/_R): 
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In order to give a transparent derivation of the solution to fl30l) and the corresponding value 
of the classical action, we first consider the two independent solutions to the homogeneous 
version of ([30a]): 



f{t) = V + sin ^Vl + aarctan— ^ , 
h{t) = VR^ + cos (^VT + a arctan— . 



(32a) 
(32b) 



A useful object to consider is the Dirichlet Green function of the operator df + aR"^/ (-R^ + 

{hh{t^)-hrf{t^)){hh{t>)-h2f{t>)) 



G{t, t'\ti, t2 



W[f,h]{hh2-hf2] 



(33) 



satisfying 



aR^ \ 

d'f + - -^]G{t,t'\ti,t2) = 6{t-t'), G{ti,t'\ti,t2) = 0, G{t2,t'\tut2) = 

{R^ + t^) J 

(34) 

with W[f, h] = fh — hf being the Wronskian of f(t) and h(t) (independent of t), t2 > ti, 
= min(t,t'), t> = max(t,t') and /i = f(ti), hi = h(ti), etc. 

With the Green function given by ( j33l) . and b{t) denoting the right hand side of equation 
f l30ap . one can write down the solution to (1301) as 



Xciit\Xi,ti; X2,t2) — —Xidt'G{t,t'\ti^t2] 



+x2dt,G{t,t'\ti,t2: 

+ J dt' G{t,t'\ti,t2)b{t') 



t'=ti 

t2 



t'=t2 



(35) 



Given the above formulas, the classical action can be written in a relatively general form 
that will turn out to be useful later: 



h2fi - f2hi 



/l/l2 - hif2 

h2f{t) - f2h{t) 



/l/l2 - hif2 
/l/l2 - hif2 



xl - 


\ W[f,h] 1 




Jih2 - hif2_ 



t2 



9 - 
- -— 

- £ dt b{t) ^ ^^^^ _ ^^^^ 

Jti Jti 

.2 



Xl - A;_ dt bit) 

'h V /i^2 - f2hl 



hh{t) - hf{t) 



W[f,h]{fih2-hj2) 



^(t2-ti)--ln. 



2 Vi^^Tt? 



X1X2 (36a) 

X2 (36b) 
(36c) 
(36d) 
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Of course, one can also solve the equations of motion fl30|) explicitly: 



^R^ + sin2Ai2 



V , - --.-^t2 , ^ VR^ + t'^ sin2Ait 

^ — ^1 , _ ^ „ : — — — h A2 



^/WT^ sin2Ai2 

sin2A(2 sin2Ait 



(3keVR^ + 



A12 

At2 



Rk^l + a) 
1 — 



1 - 



sin2Ai2 sin2A 
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y/Y+~a f , t2 , ti 
arctan arctan — 



R 



yjl + a [ t2 . t 

I QTT'rQ'n. 

2 



( arctan— — arctan— 
\ R R 

( 



t ti 

arctan arctan — 

\ R R 



(37a) 
(37b) 
(37c) 
(37d) 



The classical action can now be evaluated either by brute force using the explicit classical 
solution (1371) . or, with less work, from (l36l) : 



S4Xi,ti\X2,t2] = - 

+ k. 



u 



R^l + a 

27- COt2Ai2 



2 (i?2 + t\) 2 (i?2 + t\ 
t2 , R^yTTa 



+ 



2 (i?2 + tl) 2 {R^ + t 
k_y/l + aR 

Pke 

VTT^^/R^ + tl 

f3ke 
VTT^^/R^ + tl 



cot2A 



12 



tanAi2 
tanAi2 



12 



X1X2 



Xo 



2 U2 



A;_ (1 + a 
,2 



,3/2 



-ftanAio - A 



R^ 



i \/R^ + t 
^(t2-ti)--ln 



EA 



e^i2 



2 /R^Tt? A;_i?v^rT^' 



(38a) 
(38b) 
(38c) 
(38d) 
(38e) 
(38f) 

(38g) 



Next we consider the "quantum-mechanical" prefactor fHOl) . With the expression for Sd 
given by fl36l) . equation fl29l) becomes: 



a^(ti,t) _ i fih-hj 
at 2hh-hj 



This leads to the solution 



A{t^,t2) = AT (i?2 + t2)-l/4 ^^2 ^ ^2^-1/4 | sin2Ai2r' 



/2 



5m, 



(39) 



(40) 
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which contains the Maslov phase 0m and a constant normahzation factor M . The Maslov 
phase is piecewise constant away from the focal points (the positions of focal point for ^2 
are functions of ti). Its value is worked out in appendix [Bl 

We can now fix the normalization of A{ti,t2) by imposing fl25]l : 



3.2.3 The i? ^ limit 

Now that we have constructed the "propagator" 0(Xi, ti|X2, ^2), which provides a basis of 
mode functions labelled by Xi, we can investigate its — > limit, which will only exist for 
special values of a and (3. 

The first non-trivial condition for the existence of an i? — limit comes from the 
prefactor ^(^1,^2)- This will vanish identically for ti < 0, ^2 > as i? is sent to (which 
would make the field operator vanish identically and manifestly destroy unitarity), unless 

a = - 1, (42) 

with an integer N , on account of the structure A/-R/sin2Ai2 in fj40|) combined with fj4Tl) . 

This behavior of the prefactor ^(^1,^2) can be naturally understood by inspecting the 
classical homogeneous solutions (1521) . For generic values of a and (3, in the i? — > limit, 
those behave as f{t) —>■ t and h{t) —>■ \t\. Since these two functions are not linearly 
independent on the negative real axis, one will not be able to specify arbitrary initial 
conditions (Xi,Vi) for the classical solution at ti < 0. Should one try to do so, in the R ^ 
limit, the classical trajectory will be kicked away to infinity for all t > 0. Correspondingly, 
all quantum wavepackets will be kicked away to infinity, and the wave function will vanish at 
all finite values of X for t > 0, as manifested by the behavior of the prefactor A{ti, ^2)- This 
problem is avoided, however, for special values of a: if a = (2X)^ — 1 (with an integer N), 
the i? — > limit of the two solutions is f{t) — > sign(t) and h{t) \t\; ii a = {2N + 1)^ — 1, 
it is f{t) t and h{t) 1. 

A further condition on a and arises from considering the R ^ limit of the classical 
action Sci[Xi, ti|X2, ^2]- The problematic terms in the action fl55]) are those with coefficients 
containing namely: 



tanAi2 - A12 ) - , „® ■ (43) 



A;_(l + «)3/2i?V k.R^/T+a 
In order to cancel the divergences, we find the following equation for /3: 

3' = (44) 

tan (7rv/rTa/2) 

7ryrTa/2 
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If a = (2A^ + 1)^ — 1, this condition would make us naively conclude that (3 = 0. However, 
a direct inspection of fH3|) reveals that it wouldn't make the divergence cancel. The only 
other option remaining is 

a = {2Nf - 1, I3 = 2N (45) 

{(3 = —2N corresponds to the same space written in different coordinates). These are the 
conditions for existence of an i? — > limit for a free scalar field dynamics on the generalized 
nuU-brane. As one would expect, the values a = 3 and (3 = 2 corresponding to the original 
nuU-brane do meet these conditions (for = 1). 



3.2.4 Mode functions in momentum basis 

To facilitate comparison with earlier work, we shall now derive momentum basis mode 
functions using the position basis mode functions (Xi, ti|X2, ^2) of section [5. 2. 2[ The ex- 
istence of an i? — >■ limit will not be affected by such conversion. To obtain the momentum 
basis mode functions, we take a Fourier transform (with respect to Xi) of the "propagator" 
(Xi, ti|X2, ^2) to convert it to an incoming plane wave basis, cancel the "free evolution" 
by multiplying with exp [i [rn? + p^) ti/(2fc_)] (we are simply using the freedom we have in 
defining the basis of mode functions), and take the ti —>■ —00 limit (which refers our mo- 
mentum labels to incoming waves in the infinite past). We shall also omit time- independent 
overall phase factors. Finally, we shall take the limit ti — > —00: 

Vk_,ke,p,m = ^ lirn^ J dXi A{ti,t2) exp {-iSd) exp (ipXi) exp ( — J (46) 

To be consistent with our notation in the beginning of this section, we shall now switch 
back to writing X~^ instead of t for all space-time quantities. We finally obtain the mode 
functions 



Vk.,ke,p,m \l |^_| + (x+)2) I sin2A I 



xexp 



-ik-X^ , . IN iX ( ^ ^ 2RNp\ 



ipke ^ ip^RN ikl, ^ iim?X^ ., ., ^■ 

-tanA - cot2A + ^, ^^^ tanA + h ik^X~ + ike, © 



A;_ k^ 2k_NR 2k. 

(47) 

where A = Xarctan(X+/i?) + 'kN/2. 
For i? — > we obtain: 



1 /ivr \ 

Vk_,ke,p,m =-===exp (^-sign(fc_)(2Ar - 1)^(X+) j 



X exp 



-tpX sign(X+) + + + ik^X + ikeQ 



(4J 
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where 6{t) denotes the Heaviside step function. 

In order to compare our mode functions P7|) with those derived in [8], we write out the 
mode functions for = 1 exphcitly: 

,,(iv=i) 1 , r -ik-R^X'' iX / R + 



'k.,ke,P,m ^y\k^C'^^''l2X+{R^ + {x+y) ^i?2 + (x+)2 V^®x++^ x+ 

tpkeR ( R^-{X+f \ ikl irn^X^ ^ % n 

-kJi^-kiy J-2fc:X^ + ^fc^ + ^^-'' +^^«® 

(49) 

The general expression for the Maslov phase (^m is given in appendix [Bl For = 1, there 
is only one focal point at X^ = 0, and the Maslov phase becomes: 

0M = exp(^|^sign(fc_)^(X+)^. (50) 

Written in our interpolating coordinates from f|T6|) . the mode functions of [8] can be 
re-expressed as 



=exp 



ip+R^X^ %X fR^ + iX+f X- 

J ;ttt: n- 



./lX+ "L2X+(i?2 + (X+)2) ^R2 + (x+)2 V RX+ R 

iX+fn~jV iJ^ im'^X+ 

H ip^ X + mB 



(51) 



2p+ \ R J 2p+X+ 2p 



with J, n being the labels used in |8]. 

The equality of the two expressions (up to normalization conventions) is established by 
recognizing the following identifications: 

k_ = —p~^, 

P = (52) 
k@ = n. 

Our results thus agree with those of [8] for non-zero values of R in the particular case of 
the null-brane (a = 3, /3 = 2). Now we are ready to examine the i? — >■ limit (where our 
choice of coordinates will reveal a peculiar reflection property at the singularity) for the 
original as well as the generalized null-brane. 

3.3 Qualitative properties of the singular limit 

To recapitulate, we have examined the dynamics of a free scalar field on the following 
3-parameter (a, (3, R) family of backgrounds ("generalized null-brane"): 

ds' = -2dXUX~ + ^Tf~^h dX^r + -,M£^dXUQ 

(i?2 + (X+)2)2^ ' ^R2 + (X+)2 (53) 

+ {R^ + (x+)2) de^ + dx^ 
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As R goes to 0, all of these geometries (irrespectively of the values of a and /5) reduce (away 
from = 0) to the parabolic orbifold times a line: 

ds"^ = -2dX+dX- + (X+)2c/e^ + dX^ (54) 

What we have found is that the i? — > limit of the scalar field mode functions exists only 
when a = {2N)'^ — 1 and (3 = 2N, with N being an integer. 

In terms of the limiting expression for the mode functions fH8|) . we find few surprises. 
The result is essentially independent of the values of a and (3 (for those values for which the 
limit exists) and bears a close resemblance to the mode functions obtained in [8] . However, 
the minor discrepancy between our results and those of [8] deserves some clarification. 

For any a and (3, the metric of the generalized null-brane converges to the metric of the 
parabolic orbifold times a line (15^ . which is formally the same as the metric of the parabolic 
orbifold written in the y-coordinates of [8j. Given only the R = expressions, it may 
therefore be tempting to identify {y~^, y~ ,y, u) {X'^ , X' , 9, X). With this identification, 
however, the mode functions are not exactly the same, even for the case of the original null- 
brane {a = 3, (3 = 2). The difference between the two is the factor of sign(X"'") in front of the 
ipX term in the exponential of fj48|) . It is important to realize that the difference between 
the two sets of mode function does not represent any dynamical distinction. Rather, it is 
explained by the difference in the choice of coordinates. 

To construct the parabolic orbifold flM|) as an i? — > limit of the null-brane, the authors 
of [8] employ their singular y-coordinates (this coordinate system fails at X~^ = even for 
smooth spaces at non-zero R). As a result, they obtain mode functions withouiEl the 
aforementioned factor sign(X"'"). 

On the other hand, we construct the parabolic orbifold metric 0541) and the correspond- 
ing coordinates as an i? — limit of smooth coordinate systems parametrizing smooth 
geometries. In this case, the factor of sign(X''') is present and its effect is that the position 
and velocity in the X-direction for all particles are reflected as they pass through X~^ = 0. 

Even though the two sets of mode functions are essentially equivalent (and only differ 
by a coordinate choice), one may think of our present parametization as being more ac- 
curate. Indeed, it is very natural to demand that, since the singular space is constructed 
as an ^ limit of smooth resolved geometries, the coordinates on the singular space 
should be constructed as an i? — limit of smooth coordinate systems on the smooth 
resolved geometries (even though, with our present theoretical understanding of space-time 
singularities, it is not possible to give a systematic justification to this treatment of coordi- 
nate systems). Note that the flip of the X-direction for positive X^ cannot be undone by 
a smooth coordinate transformation, so it will always be present if the parabolic orbifold 
metric flMl) is constructed as a limit of a smooth coordinate frame on the (generalized) 
null-brane. 

One finds considerably more surprises if one contemplates the properties of those geo- 
metrical resolutions for which the singular limit of the scalar field dynamics exists (rather 

^Incidentally, the same mode functions are obtained by applying the non-geometrical "minimal subtrac- 
tion" prescription of [13j directly to a free scalar field on the parabolic orbifold, without any recourse to 
the null-brane or its generalizations. See appendix [Al 
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than merely examining the hmiting expressions for the mode functions). 

Firstly, looking at the null-brane example of [8], one could get the impression that the 
singular limit exists because the curvature is identically zero for any finite R (in a way, one 
can say that the singularity is never really "there"). However, this is not the case. The 
non- vanishing components of the Riemann tensor for our generalized null-bran^ geometries 
are 

_ R\4a-3P^) _ R^P^-4) 

- 4(i?2 + (X+)T' ~ 4(i?2 + (x+)2)- ^^^^ 

So, generically, the curvature will blow up around X"*" = as i? is sent to 0, even for 
those values of a and /3, for which the singular limit of the scalar field dynamics exists 
(a = {2NY — 1, /? = 2N). Note that the Ricci scalar vanishes, so our results do not depend 
on the choice of the Ricci scalar coupling of the scalar field (minimal, conformal, etc., [25]). 
Furthermore, one could examine the Weyl tensor: 

- 2(i?2 + (x+)2)2' - - 2{R^ + (X+)2) 

and notice that it actually does vanish for all those cases when the limit exists (within the 
particular family of geometries we have been considering). However, there are many values 
of a and (3 for which the Weyl tensor fl56l) vanishes, yet no i? limit of the scalar field 
dynamics exists. For that reason, conformal flatness is not likely to constitute an important 
part in possible explanations for the existence of the singular limit. 

Perhaps the most puzzling feature of our results is the very fact that the limit appears to 
exist for a discrete subset of the possible parameter values within our family of geometries. 
One could think of this as being an artifact of choosing our particular slice in the space of all 
possible geometries (this, however, would obviously require a fairly delicate coincidence). If, 
on the other hand, the feature is generic, it would point to an interesting sort of discreteness 
inherent to the dynamics in the near-singular region. This question would certainly deserve 
further investigation, even though that would require mathematical machinery going beyond 
what has been employed in our present considerations. 



4 The light-like reflector plane 

Our generalized null-brane considerations of the previous section suggest a very natural 
simplification. Namely, instead of 0171) we shall consider the following metric: 

ds^ = ^ (dX+Y - 2dX+dX- + dX"". (57) 

This space-time can be classified as a pp-wave geometry (see, for example, [26]), and, in 
particular, it can be extended to a 10-dimensional background satisfying Einstein's equation 

^These will obviously vanish for the values of a and /3 corresponding to the original null-brane space-time: 
a = 3, /3 = 2. 
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by inclusion of the appropriate antisymmetric tensor field and dilaton. What makes this 
particular pp-wave interesting is that it develops a strong singularity when R is sent to (for 
example, the singularity is much more dangerous than the "weak pp-wave singularities" of 
[27]). Furthermore, the R limit of the wave equation in this geometry can be explicitly 
analyzed. 

We shall not present a detailed derivation of the mode functions for (!57|) . but simply 
notice that the wave equation in this background is formally analogous to that on the gener- 
alized nuU-brane (!20l) . with f] = 0, ke = and the first term on the right hand side omitted 
(this term comes from the determinant of the generalized null-brane metric). The expression 
for the position basis mode functions 0(Xi, ti|X2, ^2) = -^(^i, ^2) exp[— iS'c/(Xi, ti|X2, ^2)] 
then follows from (|38| HOl liTj) : 



Scl [Xi, ti\X2, t2 



k^ 

k^ 



RVT 



a 



2{R^ + ti) 2(i?2 + tf) 
t2 RVTTa 



2(i?2 + t2) 2{R^ + q) 
k^^/l + aR 



A{ti,t2 



^/R^ + tl^/R^ + tl sin2Ai2 
2n 



cot2Ai2 
cot2Ai2 
X1X2 - 



XI 



m 



2k_ 



■(t2-ti) 



R^ + tlJR^ + msm2 A 
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(58a) 
(58b) 
(58c) 

(59) 



\RVl + a 

(with A12 = y/l + a (arctan(t2/-R) — arctan(ti/i?)) /2 and (pM being the appropriate Maslov 
phase). 

The R —>■ limit of (piXi, ti\X2, ^2) exists if a = K"^ — 1 (with K being an integer) and 
equals 



0(Xi,ti|X2,t2) 



1 



27r to - U 



X exp 



^ (X2 - (sign(t0sign(t2))^+^Xi)' + ^(^2 - h 



(60) 



If K is odd, the above expression merely represents free motion on Minkowski space. 
To verify this statement, one can simply check that 0(Xi, ti|X2, ^2) solves the Minkowski 
space wave equation written in light cone coordinates: 



- tc 



2k^ 



+ 



m 



2k^ 



(61) 



Despite the strength of the singularity in the i? — > limit, the free scalar field dynamics 
actually becomes identical to that on a completely flat space. 

If K is even, the motion is still free for all positive and all negative t (read X^). However, 
as the particles pass through X~^ = 0, their positions and velocities in the X-direction are 
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reflected (note the (sign(ti)sign(t2)) structure in fl60l) ). This reflection is similar to the 
one happening for the case of the generahzed nuU-brane, but it occurs on a simpler space- 
time geometry. The mode functions corresponding to fl60|) can be easily derived by means 
of a Fourier transformation, analogously to section I3.2.4[ 

Because of the property we have just described, we call the space-time flFTI) with a = 
{2NY — 1, where N is an integer, the light-like reflector plane. It is an extremely simple 
family of pp-wave geometries developing a strong singularity at = when R is sent 
to 0. Furthermore, at least for free propagation in this background, the singular limit is 
manifestly well-deflned, and includes a curious light-like object reflecting the positions and 
velocities of all particles passing through it. Given its simplicity, the background presented 
in this section may be worth studying as a toy model for light-like singularities in both 
perturbative string theory and matrix theory contexts. 

5 Conclusions 

In this paper, we have noted the relevance of Hamiltonians with singular time dependences 
involving multiple operator structures for the problem of geometrical resolution of singular 
space-times. We have given a general review of the quantum dynamics corresponding to 
this type of Hamiltonians, and emphasized some important simpliflcations, which can occur 
if the Hamiltonian possesses a flnite dimensional dynamical group. 

Turning further to a particular case of interest, we considered a 2-parameter general- 
ization of the nuU-brane space-time, and addressed the question of the singular limit of 
a free scalar fleld dynamics on this background. Surprisingly, this limit happened to ex- 
ist for a discrete subset of the possible values of the two parameters. The limiting mode 
functions are closely related to those previously obtained for the nuU-brane [8]. We have 
opted for an accurate coordinatization of the singular limit of our spaces, based on taking 
a limit of smooth coordinate systems on the smooth geometrical resolutions. In contrast to 
those employed in our coordinate system reveals a peculiar "reflection" property of the 
generalized (as well as the original) null-brane space-times. 

Our analysis of the generalized null-brane suggested a natural simpliflcation, and we 
termed this simplifled geometry the light-like reflector plane. This space-time is a fairly 
simple pp-wave, whose singular limit is a light-like plane reflecting the positions and veloc- 
ities of all particles as they pass through it. Due to the simphcity of its singular structure, 
the light-like reflector plane may turn out to be an interesting toy model for studying 
light-like space-time singularities in perturative string theory and matrix theory. 
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A Minimal subtraction for the parabolic orbifold 

It is instructive to compare our derivations of section 3 (based on geometrical regularizations 
of the singularity) with what we would have obtained by applying the "minimal subtraction" 
prescription of [13] directly to the parabolic orbifold. 

We write the action for a massive scalar field in the parabolic orbifold using the metric 



S = J dydy+dy-dz\y^\ - - • (62) 

We decompose into Fourier modes along y~ and y (with the condition 0^*^^ = (p^^^ky)'- 

j dl(pi,kyew{i^y' + ikyV)- (63) 

ky 

Now the action can be rewritten as 

S = Y.j^y^dl \y+\ \il [ct^i^kyd+ct^lky - - ( + ^i^kyCt^lk] (64) 

ky 

The equations of motion are 

^^Id^^iMy + ^ + + hky = 0. (65) 

One can deal with the constraints due to the first order nature of the light-cone formalism 
by choosing (pi^^y = \/ily~^(f>i,ky as the canonical coordinate and iii^ky = \f^hj^4'l,,ky its 
conjugate momentum. We obtain the Hamiltonian 

H = Y.jdy^dl^^ ( + ^^^^yk^r (66) 

ky 

We now apply the "minimal subtraction" scheme of [13] to the singular time {y^) depen- 
dence in (|66|) : 

1 {y^f - 



W {{y^f 



e 



2-] 



(67) 



(One could in principle add a (resolved) 5-function with an arbitrary coefficient on the right 
hand side, but we shall not make use of this freedom for the sake of brevity.) The solution 
for (j) reads: 
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We can return to the original (pi^ky and write the scalar field mode functions as: 

1 / I/+ \ 

'il^ik m?{v^ -iV) exp I tV^ -\ — tt:^ n + ikojV + 'ilv^ 1 • (69) 

ri,/c„m vy ,y ,yj ,/2ily+ \ 2ir 2il {y+Y + ^ ) ^ ' 

In order to compare with the mode functions of [8] , we identify ky = J and p"^ = —p^ = 
—I. Furthermore, the mode functions of [8] are derived for the parabolic orbifold times a 
line, whereas the mode functions (jniD refer to the parabolic orbifold proper. To compensate 
for this difference, one should set the momentum along the extra line in [S] to zero (which 
amounts to imposing n = J in the notation of that paper). Thereafter, the two sets of 
mode functions agree. 

Note that the agreement is largely accidental. The mode function of [8] are derived via 
a geometrical regularization, but they are written in the singular y-coordinates (see further 
discussion in section 3.3). The minimal subtraction mode function are obtained through a 
regularization procedure that does not admit a geometrical interpretation. 



B Focusing properties of the wave equation and the Maslov phase 

In evolution of classical dynamical systems, it often happens that all classical trajectories 
starting at {Xi,ti) reach the same point X*{Xi) at the moment t*(ti) (irrespectively of 
their initial velocity Vi). Under such circumstances, one speaks of t*{ti) as being a classical 
focal point of the evolution. 

If t*{ti) is such a focal point, the classical action Sci[Xi,ti\X2,t*(ti)] will diverge unless 
X2 = X*{Xi,ti) (for, did it not, there would have been classical trajectories connecting 
{Xi,ti) and {X2,t*{ti)) for X2 7^ X*{Xi,ti), in contradiction with the definition of a focal 
point). As we already remarked in section [3. 2. 2[ if one pursues a semiclassical construction 
of the quantum-mechanical mode functions, such singular behavior of the classical action 
introduces formal complications in equation fl29|) . A general recipe for handling this type 
of singular features can be given [20l[2T]. However, in our present context, it will be more 
practical to analyze the relevant solution for our particular form of the classical action. 

To this end, we shall rewrite fl5^ in the following form: 



Sci[Xi,ti\X,t] = — 



fih - hj 



{X-X*iX,,t,,t)f + ■■■ (70) 



fih-hf_ 

^~ ^ \n[f,h~ hJ]iX - X*{X,,t,,t)f + (71) 



2 dt 

where the dots represent contributions non-singular ai t = t* [ti) . A focal point X*{Xi,ti, t*{ti 
is reached whenever 

hh{n - hfin = fihMt*) - hit,)fit*) = 0. (72) 

At the same time, equation fl39|) for the prefactor A{ti,t) can be solved on the left and on 
the right of the focal point t*(ti) (even though constructing a solution at t = t*{ti) naively 
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would be problematic on account of the singularity on the right hand side of fl39|) ): 

Aih,t) = M> \hh{t) - hjit)\-'/' (t > r(ti)) 

(with Af^ and A/'> being complex constants). 

Armed with these relations, we can examine the behavior of the entire wave function 
(j){Xi,ti\X, t) = A{ti, t)exp {—iSci [Xi, ti\X, t]) in the vicinity of a focal point t*(ti). Gener- 
ically assuming that fih{t) — hif{t) has a simple zero at the focal point, fih{t) — hifit) oc 
t — t*{ti), and keeping in mind that 



lim \ - — ^ exp ( — — sign(A) ) exp (iAx^) = 6(x), (74) 

A^oo V TT V 4 / 



we conclude that 



with 



lim (j){Xi,ti\X,t) = A<6(X - X*(Xi,ti)), 
lim 0(Xi,ti|X,t) = A>6{X -X*{Xi,ti)), 



A<: A/'< exp (ivr sign(A;_)/4) 



(75) 



(76) 



Ay A/'>exp (— ZTT sign(A;_)/4) 

If we further demand that the limits in (175]) should be the same (this automatically ensures 
that any convolution of (j){Xi,ti\X,t) with a smooth wave packet is continuous across the 
focal point), we conclude that 

^ = exp ^ysign(A;_)^ . (77) 

When there are many focal points t|(ti), each of them will give a contribution, and the 
resulting wave function can be written as 

0(Xi, h |X, t) = Af<j)M \fih{t) - hj{t)r'/' exp {-tS,i [Xi, h\X, t]) (78) 

with a constant normalization factor Af and the Maslov phase (pM of the form 



exppjsign(fc_)5^^^(t-t*) (79) 



{9{t) being the Heaviside step function). 

We now turn to our specific case for which the relevant term in the classical action ([38 
near a focal point is given by: 



S4X.MX,t]^k_ -^^^0I±cot2A,. + (80) 
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The focal points are the poles of cot2Ait: 



t* = 



ti + Rtan{'ITe/^/T + a) 



(81) 



1 - tan(7r£/Vl + «) ti/R ' 



The value of a determines the number of focal points. We will restrict our attention to the 
case a = (2A^)^ — 1 relevant for our present investigation, and we obtain £ G {1, . . . , 2A^ — 1}, 
i.e., 2N — 1 focal points. 

In the R —y limit, all the 2N — 1 focal points are squeezed into t = 0. This means 
that there will be one large phase jump from t < to t > 0. We thus obtain the following 
expression for the Maslov phase in the R ^ limit: 



References 

[1] B. Craps, S. Sethi and E. P. Verlinde, "A matrix big bang," JHEP 0510 (2005) 005 
[ arXiv:hep-th /0506180] ; B. Craps, A. Rajaraman and S. Sethi, "Effective dynamics 
of the matrix big bang," Phys. Rev. D 73 (2006) 106005 [IarXiv :hep-th/0601062] . 

[2] T. Hertog and G. T. Horowitz, "Towards a big crunch dual," JHEP 0407 (2004) 
073 [arXiv:hep -th/0406134 ]; "Holographic description of AdS cosmologies," JHEP 
0504 (2005) 005 QarXiv : hep-th/050307l| ; N. Turok, B. Craps and T. Hertog, ; "From 
Big Crunch to Big Bang with AdS/CFT," larXiv: 0711 . 18241 [hep-th] ; B. Craps, 
T. Hertog and N. Turok, "Quantum Resolution of Cosmological Singularities using 
AdS/CFT." arXiv: 0712 .4180 [hep-th]. 

[3] D. Robbins and S. Sethi, "A matrix model for the nuU-brane," JHEP 0602 (2006) 052 
jarXiv:hep-th/0509204|; E. J. Martinec, D. Robbins and S. Sethi, "Toward the end of 
time," JHEP 0608 (2006) 025 [arXiv : hep-th/0603104j . 

[4] J. Figueroa-O'Farrill and J. Simon, "Generalized supersymmetric fluxbranes," JHEP 
0112 (2001) Oil [a rXiv:hep-th/0110170j . 

[5] G. T. Horowitz and A. R. Steif, "Singular string solutions with nonsingular initial 
data," Phys. Lett. B 258 (1991) 91. 

[6] M. Li, "A class of cosmological matrix models," Phys. Lett. B 626 (2005) 202 
[ jarXiv : hep-th/0506260| ] . 

[7] H. Liu, G. W. Moore and N. Seiberg, "Strings in a time-dependent orbifold," JHEP 
0206 (2002) 045 QarXiv : hep^th/0204168P . 

[8] H. Liu, G. W. Moore and N. Seiberg, "Strings in time- dependent orbifolds," JHEP 
0210 (2002) 031 [ .arXiv:hep-th/0206182j . 




(82) 



23 



[9] A. Lawrence, "On the instability of 3D null singularities," JHEP 0211 (2002) 019 
[ arXiv : hep-th/0205288 ] . 

[10] G. T. Horowitz and J. Polcliinski, "Instability of spacelike and null orbifold singulari- 
ties," Phys. Rev. D 66 (2002) 103512 [ arXiv : hep-th /0206228 ] . 

[11] M. Berkooz, B. Craps, D. Kutasov and G. Rajesh, "Comments on cosmological singu- 
larities in string theory," JHEP 0303 (2003) 031 [arXiv : hep-th/0212215 ] . 

[12] J. McGreevy and E. Silverstein, "The tachyon at the end of the universe," JHEP 0508, 
090 (2005) [ arXiv : hep-th/0506130i l . 

[13] B. Craps and O. Evnin, "Quantum evolution across singularities," larXiv : 0706 . 0824] 
[hep-th] . 

[14] W. Magnus, "On the exponential solution of differential equations for a linear opera- 
tor," Commun. Pure Appl. Math. 7 (1954) 649. 

[15] B. C. Hall, Lie Groups, Lie Algebras, and Representations, Springer- Verlag: New York, 
2003. 

[16] W.-M. Zhang, D. H. Feng, and R. Gilmore, "Coherent states: theory and some appli- 
cations," Rev. Mod. Phys. 62 (1990) 867. 

[17] I. A. Malkin, V. I. Man'ko. Dinamicheskie simmetrii i kogerentnye sostojanija, Nauka 
(1979). 

[18] H. P. Yuen, "Two-photon coherent states of the radiation field," Phys. Rev. A 13 2226. 
[19] V. V. Dodonov, V. I. Man'ko, and L. Rosa, "Quantum singular oscillator as a model 



of a two-ion trap", Phys. Rev. A 57 (1998) 2851 |arXiv : quant-ph/9712028P . 

[20] H. Romer, Theoretical Optics, Wiley- VCH (2004). 

[21] V. P. Maslov, Asimptoticheskie metody i teoriya vozmushchenij, Nauka (1988). 

[22] A. D. Polyanin, V. F. Zaitsev, Handbook of exact solutions for ordinary differential 
equations. Chapman and Hall/CRC (2003). 

[23] M. Fabinger and J. McGreevy, "On smooth time-dependent orbifolds and null singu- 
larities," JHEP 0306 (2003) 042 | ]arXiv:hep^W 0206196P . 

[24] D. Robbins and S. Sethi, unpublished. 

[25] N. D. Birrell and P. C. W. Davies, Quantum Fields In Curved Space, 
Cambridge University Press (1982) 



24 



M. Blau, M. O'Loughlin, G. Papadopoulos and A. A. Tseytlin, "Solvable models 
of strings in homogeneous plane wave backgrounds," Nucl. Phys. B 673 (2003) 57 
qarXiv : hep-th /0304198 ] . 

J. R. David, "Plane waves with weak singularities," JHEP 0311 (2003) 064 
[arXiv : hep-th/0303013p . 



25 



